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Abstract. - We have investigated the role of the quantum size effects in the evaluation of the force 
caused by electromagnetic vacuum fluctuations between ultra-thin films, using the dielectric tensor 
derived from the particle in a box model. Comparison with the results obtained by adopting a 
continuum dielectric model shows that, for film thicknesses of 1 4- 10 nm, the electron confinement 
causes changes in the force intensity with respect to the isotropic plasma model which range 
from 40% to few percent depending upon the film electron density and the film separation. The 
calculated force shows quantum size oscillations, which can be significant for film separation 
distances of several nanometers. The role of electron confinement in reducing the large distance 
Casimir force is discussed. 



Introduction. — Quantum mechanical forces, namely 
van der Waals and Casimir forces, between nano objects 
are expected to be important in the operation of micro 
and nano systems [1-3] . These forces depend significantly 
upon the optical properties of the interacting objects: for 
interacting films, variations in the experimental optical 
parameters, caused by film morphology, can determine a 
change in the force of the order of 10% [4, 5]. Since the 
basic work from Lifshitz and coworkers [6, 7] theoretical 
studies have been focused mainly on the determination of 
the forces between slabs, including semi-infinite slabs, on 
the basis of a continuum description of the material di- 
electric function [8-11]. This is a simplified description 
of a film, that neglects the modifications in the electronic 
structure related to the boundaries. It is expected to hold 
when the size of the film is large, so that the surfaces play 
a minor role in determining the dielectric response. For 
metallic films the calculated electronic distributions devi- 
ate significantly from the bulk behaviour when the size of 
the film is less than ten times the Fermi wavelength [12,13]. 
If the size of the film is of the order of few nanometers the 
continuum model does not provide an accurate descrip- 
tion of the film properties and boundary effects cannot 
be neglected. Such effects arise as a consequence of the 
discretization of the energy bands due to the confinement 
potential, which produces the quantization of the electron 



energy levels in various sub-bands and affects the optical 
and electrical properties [14-21]. 

The quantum model. — The characteristics features 
induced by the quantization can be described by a model 
in which independent electrons of mass m are confined 
within a distance d in the direction of the surface normal 
(particle in a box model) [22,23]. Assuming a jellium 
model and perfect planar surfaces the eigenvalue spectrum 
is simply given by: 
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i.e. described by a continuous quantum number fcii , giving 
the modulus of the parallel wavevector, and a discrete sub- 
band index n coming from the quantization of the perpen- 
dicular wavevector k±. The corresponding wavefunctions 
are: 



(2) 



where V is the volume of the quantum well given by 
the product of the well surface A and the well thickness 
d. In this simple model the electrons behave as a two- 
dimensional gas along the x and y space directions and 
as standing waves in the z direction, with nodes on the 
boundaries. As a first approximation one can assume the 
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size d of the quantum well to be the same as the size D 
of the ion distribution of the film. This approximation is 
too crude since it does not allow the spilling of the elec- 
tron density past the film boundaries given by the positive 
charge distribution, thus leading to a depletion of negative 
charge near the surfaces. One way to eliminate this incon- 
venience without introducing softer boundaries is to allow 
the electron charge to be distributed on a larger size than 
the positive charge. It is common practice to discuss the 
effects of the confining potential in deposited films using 
the phase accumulation model [24] , according to which the 
condition that is satisfied by a quantum well state is 
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where </>a and 4>b are the phases of the eigenfunctions 
accumulated at the two film interfaces and n are integer 
numbers. For a free standing film, limited by two vacuum- 
film interfaces, 4>a = 4>b- The case of an infinite quantum 
well is recovered by simply imposing 4>a — —n- For softer 
profiles of the confinement potential the phase shifts are 
expected to increase and to depend upon the energy. One 
can introduce a more realistic description of the film elec- 
tron states, still using the infinite potential well, by allow- 
ing the effective width d to be larger than the size of the 
ion distribution d = D + 2 A, thus obtaining the quanti- 
zation condition D + 2A„ = mr/kj_, where A„ gives the 
shift in the potential well width that allows to reproduce 
the charge spilling out of the n-th state. This introduces 
an energy dependence of the shift that in principle could 
be obtained by fitting the quantum well energies to exper- 
imental data or to the energy level distribution resulting 
from a first principle calculations. Previous works [25,26] 
have shown that the energy spectrum is very sensitive to 
the position of the barriers and relatively insensitive to 
the barriers height. For the present preliminary study we 
take an averaged d value obtained by simply imposing the 
film to be neutral with an electron charge of size d and 
a positive charge of size D. This leads to the following 
expression for the effective film thickness: 
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where mp = kpd/ir, hp is the Fermi wavevector that 
is related to the ion electron density No by the relation 
kp = 3tv 2 Nq, and mo is the integer part of tuf- Notice 
that, once D and the positive charge density No are given, 
the d value can be determined unambiguously 1 . For large 
D, when mo ~ mp , one obtains d = D + 37r/4fcp, which 
shows that the spilling out of the charge in this limit is 



1 ~For the calculation of the film charge density see for example 
reference [23]. 
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Fig. 1: Static value of the zz component of the dielectric tensor 
(continuous line) and ratio between electron densities with and 
without quantum confinement (dashed line) as a function of 
dkF/it for a film with Q p — 10 16 rad/sec. In both cases kinks 
occur, with discontinuities in the derivative, at kp/n multiples. 



proportional to the Fermi wavelength. 
With respect to the classical description of thin film opti- 
cal and transport properties [27,28], where surface effects 
are expressed as boundary conditions on the electron dis- 
tribution function in terms of the mean free path and the 
fraction of the electrons scattered specularly at the sur- 
face, this theory has the surface effects incorporated as 
boundary conditions in the one-electron hamiltonian, the 
only parameters being the film size and the Fermi energy. 

The model dielectric tensor. — The derivation of 
the dielectric function of the film described by the previous 
model, leads to a description of the optical properties that 
accounts for sub-bands transitions [22]. The finite exten- 
sion of the film along the z direction implies the anisotropy 
of its dielectric response with the diagonal components: 
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here a — x,y,z labels the cartesian component of the 
tensor, p a indicates the component of the electron linear 
momentum, lu p — Cl p N/No is t ne plasma frequency of the 
quantized electron gas (Q p — \J 'Aire 2 N / 'm is the free elec- 
tron plasma frequency) and f(Ek,. n ) is the occupation 
factor of the (k|| , n) state. The off-diagonal component 
are equal to zero. This expression differs from the plasma 
model dielectric function adopted in previous studies in 



that: (i) it has a tensor character with e x 



(ii) the plasma frequency lo p depends upon the film den- 
sity N, which changes as a function of the film thickness, 

(iii) through the double sum in the second member it ac- 
counts for transitions between lateral sub-bands, whose 
probability amplitude is expressed by the momentum ma- 
trix element between the one electron wavefunctions (2). 
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Fig. 2: Two identical interacting thin films, definition of £, D 
and the electronic spill-out A. 



It can be easily shown that these transitions do not af- 
fect the lateral components of the dielectric tensor, which 
are given by the simple expression of the plasma dielectric 
function 



,(w) = e yy (w) 
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because the momentum matrix element for x and y com- 
ponent vanishes. They modify e zz (o>) whose low frequency 
behaviour gives a finite dielectric constant: 
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here ao is the Bohr radius, and 
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The metallic character, with vanishing minimum excita- 
tion energy and divergent e zz {uj) at low frequency, is re- 
covered for d going to infinity. The plot of e zz (0), given in 
fig. 1 , shows an oscillatory behaviour as a function of d with 
periodicity given by it /hp superimposed over a regularly 
increasing curve. These oscillations arise from the peri- 
odic crossing of the Fermi energy by the sub-bands with 
increasing thickness. At each crossing a new sub-bands 
gets filled and optical transitions having such sub-band 
as initial state become possible. The occurrence of quan- 
tum size oscillations has been pointed out by several au- 
thors specially with reference to the electron density, the 
total electronic energy and the film electrical conductiv- 
ity [23,29,30]. Fig 1 shows that they are distinct features 
of the dielectric response too. 

The force. — The expression of the force per unit 
area F at T = 0° K in the configuration illustrated in the 
inset of fig. 2 can be obtained by extending previous results 
relative to isotropic films [8,31] to the case of films with 
anisotropic dielectric tensor. One obtains: 
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This equation, which to the authors khowledge has never 
been used previously for vacuum fluctuation force calcu- 
lations, deserves a few comments, (i) it has been obtained 
by considering the zero point energies associated with the 
electromagnetic modes of films of finite thickness under 
the assumption that the dielectric permittivity is repre- 
sented by an anisotropic diagonal tensor. It differs from 
the original Lifshitz formula both because it depends upon 
the film size D (i.e. has been obtained with the electro- 
magnetic field boundary conditions appropriate to a finite 
size film and not for a semi-infinite system) and for the 
presence of the anisotropic permittivity, (ii) the force goes 
to zero as the film size D vanishes, since both Qtm and 
Qte go to zero, (iii) to calculate the force one has to 
determine the frequency dependence of the dielectric ten- 
sor, which we have obtained from equation (6) using the 
eigevalues and the wavefunctions of the quantized film. 
This approach differs from what is commonly done in 
dispersion forces calculations based on a continuum de- 
scription of the dielectric properties, where the dielectric 
function is derived from empirical expressions which do 
not bear a direct relation to the electronic band struc- 
tures of the interacting bodies [8-11]. In the following 
the force per unit area calculated for the dielectric tensor 
(6) appropriate to the film is denoted by Fq, while Fp 
indicates the force calculated using the isotropic plasma 
model i.e. e xx = e yy = e zz . To evaluate the importance of 
the size quantizations, we are interested to compare Fq, 
calculated for given film thickness D and ion density No 
at different distances £, with Fp calculated in the same 
configuration. Fig. 3 (a) displays curves of Fq calculated 
at fixed Sip and £ values as a function of the film thick- 
ness. The fi p values correspond to free carrier densities of 
heavily doped semiconductors. These systems are those 
which display higher modifications. We have considered 
distances I ranging from 10 to 50 nm. For comparison the 
curves of the corresponding Fp vales are reported. One 
can see significant differences between the two models for 
d values that are of the order of few multiple integer of the 
half of the Fermi wavelength. The curves show quantum 
size oscillations with the expected periodicity kpd = mr 
superimposed over a regularly increasing behaviour. For 
thick films the results of the isotropic plasma model are 
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Fig. 3: (a) Black lines represent Fq as a function of the film 
thickness for different film densities and separations, Q. p = 10 14 
rad/sec and £ — 10 nm (continuous line), Q p = 5 • 10 1 rad/sec 
and £ = 50 nm (dotted line), Q p = 10 14 rad/sec and £ = 50 
nm (dashed line). The gray dashed lines represent Fp for the 
same parameters, (b) relative percentual difference between 
the force with and without quantum size effects. 
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Fig. 4: (a) Black lines represent Fq as a function of the plasma 
frequency for different film thicknesses and separations, D = 1 
nm and £ = 10 nm (continuous line), D = 5 nm and £ = 10 
nm (dotted line), D = 5 nm and £ — 50 nm (dashed line). 
The gray dashed lines represent Fp for the same parameters, 
(b) relative percentual difference between the force with and 
without quantum size effects. 



recovered. The size induced modifications are better illus- 
trated by the quantity: 



Fp - Fq 
F P 



(15) 



which gives the relative variation of the force with respect 
to the isotropic plasma model. Fig 3 (b) presents plots of 
8 for the cases under consideration. It is seen that the size 
induced modifications are very large, ranging from 50% to 
10%, for films of nanometric thickness even at distances I 
of several nanometers. 

One can see similar modifications in the theoretical re- 
sults obtained by keeping the film size constant and chang- 
ing its density. This is illustrated in figs. 4 (a) and (b) 
which display force curves obtained as a function of plasma 
frequency. Notice that even for typical metallic densities 
(fip of the order of 10 15 10 16 rad/sec) the deviations 
from the plasma model are quite significant and can be of 
the order of several percents for film separation distance 
I ranging from 10 to 50 nm. In a previous study on the 
Casimir effect for metal and semiconductor slabs it has 
been pointed out that the Casimir force can be consid- 



erably reduced by decreasing the slab thickness [10]. It 
is interesting to see whether quantum confinement effects 
are important in determining the reduction at large t val- 
ues. To this end we plot in fig. 5 (a) the reduction factors 
VP = Fp/Fcas and tj q = Fq/F C as [32], where 
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is the force between perfectly reflecting mirrors at separa- 
tion distance I [33] . We consider films of nominal thickness 
D = 1 -j- 5 nm and different electron densities. One can 
see that the effects of the quantum confinement tend to 
decrease the reduction factor over a large distance range. 
The correction depends upon both the film size and the 
electron density, being larger for smaller values of these 
quantities. Fig. 5 (b) plots the calculated relative differ- 
ence 6 as a function of I for the same cases. It shows that 
relative variations of the force remain significant even at 
very large distances thus giving sizeable reduction to the 
force in the large t Casimir regime. Notice that 5 does 
not show a monotonous behaviour as a function of t: it 
decreases regularly at short distances, it remains constant 
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Fig. 5: (a) Black lines represent t/q as a function of the films 
distance £ for different densities and thicknesses: continuous 
line f2 p = 10 16 rad/sec D = 5 nm, dashed line f2 p = 10 15 
rad/sec D = 1 nm, dotted line S7 P = 10 15 rad/sec D = 5 
nm. Gray lines represent ??p as a function of I for the same 
parameters, (b) Percentual difference (15) for the three cases 
of the (a) plot. 



in a large range of £ values and, before going to zero for 
very large distances, it shows a local maximum that shifts 
at higher £ values upon decreasing the thickness and/or 
the electron density. This behaviour is quite typical and 
can be reproduced in several cases. To understand this 
behaviour we notice that S can be written as 



6 = 1-™ 
VP 



(17) 



Fig. 6 shows the typical behaviour of the numerator and 
of 1/rjp as a function of £. The first curve decreases more 
or less regularly upon increasing the film separation. The 
change in the slope arises from different inverse power be- 
haviour of the short range (dispersion force) compared to 
long range Casimir force [6-8,31]. The second curve shows 
a regular increase at short distances, it has a maximum 
and goes to zero at large £ values. Therefore the maxi- 
mum in the S corresponds to the range of £ values where 
the reduction of the Casimir force is larger compared to 
the ideal case, while the plateau is due to the combined 
effect of the decrease of r\p — t\q and the rise of the l/r)p 
curve. 




Fig. 6: Continuous line represents the percentual relative dif- 
ference S, dashed line represents rjp — riQ, dotted line represents 
the ratio 1/r/p for the case H p = 10 15 rad/sec D — 5 nm 



Conclusions. — Using the particle in a box model we 
have presented results that show under what conditions 
the electron confinement can affect the quantum mechan- 
ical force between free-standing films and we have given 
evidence of the presence of quantum size oscillations effects 
for small thickness and/or low density films. In extending 
this theory to the more realistic case of deposited films, 
one has to account for the penetration of the electron 
wavefunctions into the substrate. This requires a mod- 
ification of the quantization condition compared to the 
abrupt infinite barrier potential adopted in the present 
paper. However for semiconductor substrates where the 
confinement is essentially due to the fundamental gap, the 
modifications of the potential should not alter appreciably 
our main conclusions. For metallic substrate the extension 
is less obvious since the confinement is related to the ex- 
istence of symmetry or relative gaps of the projected bulk 
band structure, that exist for particular directions only. 
In this case the solid-film interface has to be simulated 
with a softer confinement potential [34]. 
In view of the sensitivity of our results to the film density, 
the observation of the effects we have reported could be 
made possible by some experimental technique that allows 
to modify the carrier density in overlayers (doping [35], 
creating electron-hole plasma by illumination [36] etc.). 
Also by changing the substrate one may affect the confine- 
ment potential, making it more abrupt or more rounded, 
and enhance or reduce quantum size effects. 
Beside including wavefunction penetration into the sub- 
strate, theoretical improvements to the present theory 
could come from the inclusion of (i) local field effects in the 
description of the film dielectric response [37], (ii) relax- 
ation time into the metal dielectric function following the 
prescription appropriate to finite size systems [38]. Work 
along these lines is in progress. 
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